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Abstract 

The DTS experiment is a spherical Couette flow experiment with an imposed dipolar magnetic 
field. Liquid sodium is used as a working fluid. In a series of measurement campaigns, we have 
obtained data on the mean axisymmetric velocity, the mean induced magnetic field and elec- 
tric potentials. All these quantities are coupled through the induction equation. In particular, a 
strong w-eff'ect is produced by differential rotation within the fluid shell, inducing a significant 
azimuthal magnetic field. Taking advantage of the simple spherical geometry of the experiment, 
I expand the azimuthal and meridional fields into Legendre polynomials and derive the expres- 
sions that relate all measurements to the radial functions of the velocity field for each harmonic 
degree. For small magnetic Reynolds numbers Rm the relations are linear, and the azimuthal 
and meridional equations decouple. Selecting a set of measurements for a given rotation fre- 
quency of the inner sphere (Rm ^ 9.4), I invert simultaneously the velocity and the magnetic 
^ \ data and thus reconstruct both the azimuthal and the meridional fields within the fluid shell. The 

results demonstrate the good internal consistency of the measurements, and indicate that turbu- 
^ I lent non-axisymmetric fluctuations do not contribute significantly to the axisymmetric magnetic 

induction. 

Resume 

Cartographie de I'induction magntique dans un ecoulement de Couette spherique soumis 
a un champ magnetique. L' experience DTS consiste en un ecoulement de Couette spherique 
soumis a un champ magnetique dipolaire. Le fluide utilise est du sodium liquide. Au cours 
d'une serie de campagnes de mesure, nous avons obtenu des donnees sur le champ de vitesse 
moyen axisymetrique, le champ magnetique moyen, et le potentiel electrique. Toutes ces quan- 
tites sont couplees a travers F equation d' induction. En particulier, la rotation differentielle du 
fluide produit un fort effet w qui induit un champ magnetique azimutal consequent. Profitant 
de la geometric spherique de I'experience, je developpe les champs azimutaux et meridionaux 
en polynomes de Legendre et j'obtiens les expressions qui relient toutes les mesures aux fonc- 
tions radiales du champ de vitesse pour chaque degre. Pour de petits nombres de Reynolds 
magnetiques Rm les relations sont lineaires et les equations azimutale et meridionale sont decouplees. 
Je selectionne un jeu de mesures pour une vitesse de rotation donnee de la sphere interne (Rm ^ 
9.4) et j' inverse simultanement les donnees de vitesse et magnetiques, reconstruisant ainsi a la 
fois les champs azimutaux et meridionaux dans la coquille fluide. Les resultats demontrent la 
bonne coherence des mesures et indiquent que les fluctuations turbulentes non-axisymetriques 
ne contribuent pas de fa9on significative a I'induction magnetique axisymetrique. 
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1. Introduction 



It is now well establi shed that the magnetic field of most planets and stars is generated by 



the dynamo mechanism (Larmoi 



trie magnetic i ield or most planets and stars is generated by 
1919 : Elsasser . 119461: iParkerl 119551) . in which motions within 



an electrically conducting medium amplify infinitesimally small magnetic field fluctuations up 
to a level where the Lorentz force that results is large enough to stop their amplification. This is 
possible for large enough values of the magnetic Reynolds number Rm - UL/j] (where t/ is a 

typical flow velocity, L a typical length, and is the magnetic diff usivity of the medium). 

After the success of the first fluid dynarno experiments in Riga dOailitis et aL , 2001 ; Gailitis et al 

2008b and Karlsruhe dStieghtz and Miilleii I2OOII: iMueller et all l2008b . it was felt fliat the next 
step was to obtain dynamo action in a highly turbulent free-to-adjust flow. Indeed, the flow was 
very much constrained by walls and pipes in those pioneering liquid sodium experiments, while 
the flow has much more freedom in natural dynamos, and the amplitude of turbulent fluctuations 
was less than 10% that of the mean flow, while it can be of order one or more in astrophysical 
objects. 

The challenge wa s soon addressed by seve ral teams, in Russia (Ipr ick et al 



1 ne cnaiienge wa s soon addressed py seve ral teams, in Kussia (irr icK et ai. 
dLathrop et al.L I2OOII: lO'Connell et all boOlb . and in Europe (Cardi n etalT 



20021). all using liquid sod ium as a working fluid (see recent reviews bv lVerhi 



20011) in the USA 
2002l:lMarie etal 



leetal.1 (1201 0') and 



Lathrop and Forest! (1201 Ih ). In several of these experiments, the mean flow was well character- 
ized, and was expected to yield a dynamo above a critical magnetic Reynolds number that was 
achievable. It therefore came as a surprise that none of these experiments produced a dynamo, 
although a rich vari ety of dynamo behaviours have been d i scovered in the VKS e xperiment in 



Cadarache, France (IBerhanu et al.L 120071: iMonchaux et al.L l2007t lAumaitre et al.L |2008) when 



ferromagnetic disks stir the fluid. 

It appears that turbulent fluctuations have a collective contributi on to the mean magnetic in- 



duction that counteracts that of the mean flow dSpence et al.L 12006b . That small-scale turbulent 
fluctuations can contribute to a lar ge-scale magnet i c field is not surprising. In fact, it is the 
basis of the or-effect introduced by iSteenbeck et al. (Il966l) . who showed that the interaction of 
small-scale velocity fluctuations with small-scale magnetic field fluctuations produce an electro- 
motive force proportional to a large-scale magnetic field. The success of the Karlsruhe dynamo 
experiment relies precisely on this collective effect. 

In rotating bodies such as planets and stars, it has long been recognized that the combination 
of an a^effect with an fa^ -effect provide s an ap pro priate recipe to p roduce a large-scale magnetic 
field (|ParkeiilT955b ( see ICharbonneaul hOOSh and lRieutordf(l2008l) for reviews). The w-eff'ect is 



due to the shear caused by differential rotation in the fluid. It converts poloidal magnetic field 
into toroidal magnetic field. One then invokes some kind of c-effect, due to non-axisymmetric 
motions, to convert toroidal field into poloidal field, in order to sustain the dynamo . Many mean 



field dynamo models for the Earth or the Sun are based on this mechani sm (e.g., iParken, 11955 
Barenghi and JonesL 1 199 It Brandenburg et allll99ll : ICharbonneaull2005b . 
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It is therefore of some interest to evaluate these two effects in various situations, in particular 
to understand when, where and why turbulent fluctuations favor or hinder d ynamo action. There 
has been several efforts in this direction, using either n umerical simulations ( Spence et all 2008 ; 
[Brandenburg et aL , 2010l) . or laboratory experiments ( Petrehs et al. , 2003 : Spence et al. , 2006f) 
by observing the magnetic response of the device when a weak external magnetic test field is 
applied. The VKS team in particular has obtai ned clear evidence fo r these effects (and others) 
in von Karman flows, as nicely summarized by Verhille et al. ( 2010l) . In a similar geometry, the 
Madison group has measured a dipolar component o f the induced magnetic field, which cannot 



be explained by the mean flow jSpence et al.L l200a). More recently , a dir ect measurement of 



the ff-effect has been obtained on the same set-up (iRahbarnia et al.L l2012h . using a dedicated 



probe that measures locally the electromotive force. Turbulent fluctuations can also enhance the 
magnetic diffusivity of the fluid, yielding the so-c alled ;6-effect. Impressive measurement s of this 
effect have been obtained by the group in Perm (" prick et aL , 2010; Nosk ov et al. , 2012 ) for the 
transient flow of sodium in a torus abruptly stopped. In these studies, only a global averaged 
response of the system is obtained, with no information on the spatial distribution of the a- and 
w-effects. 

Several teams have also been able to combine experimental measurements of the induced 
field with numerica l simulations of the f low backed by veloc ity measurements on a water model 



of the experiment (ISpence et al.L l2008t iRavelet et al. ■ l2012h . The present article presents one 



effort in the same direction, using the DTS magnetized spherical Couette flow experiment. The 
DTS experiment is the only experiment so far that combines rotation and a s trong applied dipola r 



magnetic field. It was built to explore the resulting magnetostrophic regime (ICardin et al.Ll2002h 



A number of interesting observations have akeady been presented , relating to the mean flow 
dNataf et all l2006l bOOsi: iBrito et all l201lb and to the fluctuations dSchmitt et al.l l2008l l2013b . 
Since the imposed magnetic field governs the dynamics of the system, it is not possible to rely 
on a water model. However, we already know that turbulent fluctuations are strongly a f fected 



by the presence of a large scale magnetic field (INataf and Gagnierd l2008t ISchmitt et al.L 120 1 3 



Figueroa et al.Ll2012l) . and it would be most interest ing to find out how this translates into the a-. 



yS- and w-effects for this kind of flow. Fortunately, iBrito et al.l (1201 Ih have shown that velocity 
profiles measured by in situ ultrasound Doppler velocimetry provide excellent constraints of the 
mean flow in DTS . Furthermore, the simple spherical geometry of this experiment makes it 
possible to use tools developed in the context of celestial dynamos. 

Thus, my aim here is to set the stage for further steps that should enable us to map the a-, 
j0- and w-effects for the magnetized spherical Couette flow, using measurements from the DTS 
experiment. The first step consists in solving the kinematic problem of the large-scale fields, 
ignoring the contribution of the small-scale fluctuations. For future reference, I write down 
the equations and the inversion procedure that apply when only the large-scale fields are taken 
into account. I then apply this procedure to a selected set of actual DTS measurements. The 
results demonstrate the good internal consistency of the measurements and show that small-scale 
fluctuations play a negligible role. 

The organization of this article is as follows: Section |2] describes the DTS experiment, the 
measurements and the way they are processed. Section [3] explains why we treat separately the 
azimuthal part and the meridional part of the fields. The forward problem of predicting the 
induced field and related quantities from the mean flow is treated in section|4]for the azimuthal 
fields, and in section |6] for the meridional fields. Section |5] presents the inversion results for 
the azimuthal fields, and section Q for the meridional fields. Conclusions and perspectives are 
discussed in section[8] 
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Figure 1 : Sketch of the DTS set-up. In this northern quadrant of a meridional cross-section, liquid sodium is enclosed 
between the inner and outer solid shells. The inner shell is made of copper and extends from radius r, (thin Hne) to r, 
(thick Hne). It contains a strong magnet, which produces an axial dipolar magnetic field, as shown by the blue field lines. 
The fluid is entrained by the rotation of the inner sphere at frequency / around the vertical axis. The outer shell is made 
of stainless steel and extends from radius (thick line) to r,, (thin line). It is at rest in the laboratory frame. 



2. Experimental 

2.1. Experimental set-up 

The DTS experiment is a spherical Couette flow experiment with an imposed dipolar mag- 
netic field. Liquid sodium is used as a working fluid. It is sketched in figure ITT] It is contained 
between an inner sphere and a concentric outer shell, from radius r* = 74 mm to r* - 210 
mm (the superscript * refers to the dimensional quantities, and is dropped after adimensional- 
ization by r* as given in [Appendix A| i. The inner sphere consists of a 15 mm-thick copper shell 
(dimensionless radius r, to r,), which encloses a permanent magnet that produces the imposed 
dipolar magnetic field Ba, whose intensity reaches 175 mT at the equator of the inner sphere. 
The stainless steel outer shell (dimensionless radius r,, to r^) is 5 mm thick. The inner sphere can 
rotate around the vertical axis (which is the axis of the dipole) at rotation rates / - InQ. up to 
30 Hz. Although the outer shell can also rotate independently around the vertical axis in DTS , 
I only cons i der he re th e case w hen the outer sphere is at rest. Additional details can be found in 



Nataf etal.1 (l2008h and lBrito et al.. (.2011.) 



2.2. Dimensionless numbers 

In this article, I concentrate on the magnetic induction equation alone, which is governed 
by a single dimensionless number: the magnetic Reynolds number Rm, which measures the 
induction over diffusion ratio. However, the actual flow depends of course on other dimensionless 
numbers that determine the relative importance of forces in the Navier-Stokes equation. Three 
independent dimensionless numbers can be constructed, and we choose the magnetic Reynolds 
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number Rm, the Elsasser number A and the magnetic Prandtl number Pm, whose expressions 
are given in table [T] together with their values in the DTS experiment for rotation rates / = 3 Hz 
and / = 30 Hz. 



/ 


Rm 


A 


Pm 


Re 


Ha 


Lu 


Hz 






v/t] 








3 


9.4 


3 10-^ 


7.4 lO-*" 


1.3 10'' 


200 


0.5 


30 


94 


3 10-3 


1.3 10^ 



Table 1: Relevant dimensionless numbers and their values for two rotation rates / of the inner sphere. 



The magnetic Prandtl number simply compares the kinetic and magnetic diffusivities, and 
is a physical property of the fluid used (here liquid sodium). The Elsasser number has been 
introduced to compare the Lorentz and the Coriolis forces. It looks small in table [T]because we 
use Bo, the intensity of the imposed magnetic field at the equator of the outer sphere, to scale the 
magnetic field, but one should keep in mind that the actual imposed magnetic field is 23 times 
larger at the equator of the inner sphere. 

The three other dimensionless numbers listed in table [T] can be deduced from the previous 
ones. The Hartmann number Ha indicates that, even at the outer boundary, magnetic forces 
dominate over viscous forces. The Lundquist number Lu compares the time-scale of Alfven 
waves to their life-time. It is needed to convert magnetic energies to the same units as kinetic 
energies. 



2.3. Data processing 

We want to investigate the interaction of the flow with the magnetic field . We thus want to 
combine measurements of the velocity field and of the induced magnetic field. iBrito et alJ (120 lib 
have recently shown that it was possible to retrieve the large scales of the mean axisymmetric ve- 
locity field in DTS by inverting criss-crossing ultrasound Doppler velocity profiles. The induced 
magnetic field is measured at several latitudes on the outer sphere. Since the toroidal magnetic 
field vanishes at the surface, it is crucial to also measure the induced magnetic field inside the 
fluid. Additional information on the interaction of the flow with the magnetic field is obtained 
by measuring the electric potential at th e surface of the ou ter sphere, and the torque applied by 
the fluid on the spinning inner sphere. iBrito et aL ( 2011 ) present these various measurements 
and display their evolution with the rotation frequency / of the inner sphere up to 30 Hz. Here 
1 concentrate on the case with / = ±3 Hz (Rm ^ 9.4), for which we have the best velocity data 
coverage. 1 give below further indications on the processing of the data I invert. 



2.3.1. ultrasound Doppler velocity profiles 

Ultrasound Doppler velocimetry provides measurements of the component of the fluid ve- 
locity along the shooting line of the ultrasound beam, from the Doppler shift of back-scattered 
acoustic energy by small heterogeneities in the fluid. In the DTS set-up, u ltrasound transduc- 
ers can be p laced in ports at 4 different latitudes (10°, -20° and +40°) ( see iNataf et alJ (l2008h 
and lBrito et al. (201 1) for more details). Purely radial shots provide radial profiles of the radial 
velocity. Shooting at some angle away from the radial direction, we get profiles that record a 
combination of meridional and azimuthal velocities. We observe that azimuthal velocities are 
one order of magnitude larger than the meridional velocities (we record a maximum azimuthal 
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Figure 2: Data coverage in tlie northern quadrant of the (s,z) plane (s: cyUndrical radius, z: vertical axis), left: azimuthal 
fields. The colored curves are the projections in the (s,z) plane of the oft'-meridian ultrasound beams used to measure 
the azimuthal velocity of the fluid. The colored dots are the positions where the azimuthal magnetic field is measured. 
The electric potential is measured in the outer shell at latitudes given by the black squares, right: meridional fields. 
The straight lines are radial ultrasound shots. The broken straight line is a meridional shot, and the curved lines are the 
projections in the (.v,z) plane of the off-meridian shots. The radial magnetic field is measured at the cyan triangles both 
inside the fluid and at the surface. The orthoradial magnetic field data are from the magenta triangles. 



velocity of 1.8 m/s for / = ±3 Hz, whereas the maximum meridional velocity reaches 0.16 
m/s). In a first approximation, the angula r velocity in the fl uid is thus simply the fluid velocity 



measured along the ultrasound beam (e.g. jBrito et al. . 200 ih 



However, we note that the shape of the profiles depends slightly upon the direction of rotation 
of the inner sphere. This is because, while the sign of the azimuthal velocity is reversed in this 
operation, the meridional flow keeps the same sign. Therefore, 1 construct the sum and the 
difference of profiles shot for two opposite rotation frequencies / and -/. 1 thus obtain refined 
azimuthal velocity profiles from the diff'erence, and meridional profiles from the sum. 

The resulting spatial data coverage in the (s,z) upper quadrant is shown by the colored lines 
in figure lZXTI for the azimuthal velocity (left) and for the meridional velocity (right) (we will see 
in section |3]that the meridional and azimuthal fields are decoupled under our approximations.) 

In this study, I focus on time-averaged properties. For a given rotation frequency /, we record 
a series of profiles over a time lapse T. For each point of the profiles, the mean velocity is the 
time-averaged value, and the standard deviation of the mean t,,,^,,,, is given by T,„ean - t/ ^JT f, 
where r is the root mean square deviation, and T f is used as an estimate of the number of 
statistically independent samples. I retain only one every 4 points along the ray because the 
original profiles were over-sampled. 

2.3.2. magnetic field inside thefiuid 

The magnetic field is measured by Hall magnetometers placed on a narrow rectangular board 
inside a cylindrical stainless steel sleeve, which penetrates the fluid radially. The sleeve can be 
placed at three different latitudes: 10°, -20° and ±40°. The azimuthal field is measured at 6 
different radii along the board, while the radial and orthoradial components are obtained at 2 
radii. We emplace only one sleeve at a time in order to limit the perturbation of the flow that it 
produces. As for the Doppler profiles, the collection of data I invert was acquired in different 
campaigns, with / ranging from 2.7 to 4.0 Hz. All data are scaled by their actual / before being 
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combined. The resulting spatial data coverage in the (s,z) upper quadrant is shown by the colored 
symbols in figure lZXTI for the azimuthal field (left) and for the meridional field (right). 

The mean azimuthal field is deduced from a time-average of the records for a given / over a 
time-lapse of typically 300 turns, and its standard deviation from the root mean square deviations 
divided by an estimate of the square root of independent samples as described in the previous 
section. The zero reference is the average value measured when the inner sphere rotates very 
slowly (in order to average over small deviations away from a perfect dipole, which amount to 
about 1% of the dipole field). 

The induced meridional magnetic field is one order of magnitude smaller than the azimuthal 
field. For / = ±3 Hz, the maximum measured meridional field is 0.6 mT, to be compared with 
5.6 mT for the azimuthal field. The amplitude of the average field is comparable to that of its 
fluctuations. It can also be contaminated by a projection of the much stronger azimuthal field 
if the alignment is not perfect. As for the Doppler profiles, I take advantage of the differing 
symmetries of the azimuthal and meridional fields and correct for misalignment by making the 
values for the latter coincide when the sign of / is reversed. The misalignments are found to 
range between 1° and 6°. I then deduce the mean meridional field from a time-average of the 
records for a given / over a time-lapse of typically 300 turns. The zero value remains difficult 
to establish because the value measured when the inner sphere is at rest varies by up to +0.2 mT 
depending on the position in which it stopped, due to the deviations of the imposed field from 
a perfect dipole. Therefore, I assume that the induced field varies linearly with / for the small 
values of / that I consider and obtain the value for / = +3 Hz by rescaling the difference between 
the average signal at / = ±3 Hz and that at / ^ 0.7 Hz. The induced meridional magnetic field 
is thus not very precisely measured in our experiment, which is reffected in its relatively large 
error bars. 



2.3.3. magnetic field at the surface 

The three components of the magnetic field are measured at several latitudes from -57° to 
57° at the surface of the outer sphere. The magnetometers are Giant Magneto Resi stive (GMR) 
devices arranged on an in house designed electronic board (see lSchmitt et alJ (120 13h for details). 
The time-averaged azimuthal magnetic field vanishes at the surface of the fluid. Therefore, the 
surface measurements only constrain the meridional magnetic field. The actual spatial data cov- 
erage in the {s,z) upper quadrant is shown by the colored symbols at the surface of the sphere in 
figure I23TT] (right). The maximum induced surface magnetic field we measure at / = 3 Hz is 
0.04 mT, which is diflicult to measure. In particular, one must take into account the effect of tem- 
perature on the signals. We correct for these, using platinum resistance thermometers installed 
on the GMR board, and calibration curves measured under controlled temperatures and imposed 
magnetic fields. The high latitude radial probes are saturated by the dipole field and not used 
in the present analysis. The zero value remains difficult to establish because we cannot use the 
value measured when the inner sphere is at rest, since it varies by up to 0.03 mT depending on 
the position in which it stopped, due to the deviations of the imposed field from a perfect dipole. 
Therefore, I assume that the induced field varies linearly with / for the small values of / that 
I consider and obtain the value for / = +3 Hz by rescaling the difference between the average 
signal at / - +3 Hz and that at / ^ 1 Hz. The values obtained from the two hemispheres and for 
the two signs of rotation are consistent and are combined to yield my final data. 
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2.3.4. electric potentials at the surface 

Electric currents are produced through the interaction of the fluid flow with the magnetic field. 
Through Ohm's law, the electric field carries some information about these currents, together 
with the induction term U x B. The time-averaged potentials do not vary with azimuth, and thus 
relate only to the azimuthal flow. I therefore include in my azimuthal inversion the differences 
between the electric potential measured at 5 different latitudes on the outer shell, as indicated 
by the squares in figure |23T| (left). The time-averaged values and their standard deviations are 
deduced from 640s-long records. The maximum electric potential difference between electrodes 
10° apart in latitude is 0.84 mV for / = +3 Hz. 



2.3.5. torque on the inner sphere 

We will see that the magnetic torque applied by the fluid on the inner sphere is directly related 
to the strength of the induced azimuthal magnetic field at its surface. We expect the magnetic 
torque to dominate over the viscous torque on the inner sphere (whereas the opposite holds on 
the outer sphere). I therefore include this data as a further constraint in the azimuthal inversion. 
The torque is recorded for all runs, as given by the electronic drive of the motor. We substract 
from the reading t he mechanical torq ue estimated for the given rotation rate from runs without 
liquid sodium (see lBrito et alJ (1201 ih for further details). The torque is equal to -1.2 + 0.1 Nm 
at / = 3 Hz, which translates into -0.29 + 0.03 in our dimensionless units (see Appendix Ai. 



3. Field decomposition 

The problem that I consider here is purely kinematic: we do not solve for the dynamics but 
we want to relate the velocity field to the magnetic field it interacts with, using the induction 
equation: 

^ = Vx(UxB) + j]AB, (1) 
of 

where rj is the uniform magnetic diffiisivity of the fluid. 

Both the velocity field U and the magnetic field B are divergence-free and can therefore be 
decomposed into their poloidal and toroidal components. I use spherical coordinates (r, 6, ip). 
Since I only consider axisymmetric mean fields, which do not depend upon the azimuth ip, the 
decomposition can be further simplified into the following expressions (e.s.. iRoberts. (.2007 )) : 

U -up + U^e^ -V xue^ + U^e^, (2) 

where u is the potential of the meridional flow and f/^ is the azimuthal velocity field. The 
magnetic field is given by: 

B - Bp + b^e^ - Bd + bp + b^e^ = V x (A,; + a)e^ + b^e^, (3) 

where A^i is the potential of the imposed dipolar field B^, a the potential of the induced meridional 
magnetic field bp, and b^p the induced azimuthal magnetic field. 

All fields are made dimensionless, using for a time scale and r* for a length scale. The 
imposed dipolar magnetic field is scaled by its intensity Bg at the equator of the outer sphere. 
The expression of Bj is thus: 
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Bd -V X AdC^ = -7 (2 cos 9er + sin eg) . (4) 

The induced magnetic fields bp and b^p are scaled by RmBg, where Rm = Q(r*)^/77 is the 
magnetic Reynolds number, since we expect the induced field to vary almost linearly with Rm 
and with Bg in the moderate magnetic Reynolds number regime that we consider 

Led by our experimental measurements, I further consider the following approximations: 
» Up and 8^ » » bp. In other words, the flow is predominantly azimuthal, and it 
induces a magnetic field that is also mostly azimuthal, and which remains small compared to the 
imposed dipolar field. 

Under all these conditions, and assumi ng steady state, the induction equation decom poses 



into two independent scalar equations (e.g.. lFearn et al.Lll988t iBarenghi and JonesLll991h : The 



induction equation for the azimuthal field b^ yields: 

sBd ■ + Aib^ = 0, (5) 

with Ai - - l/i^, and where we introduced the angular fluid velocity a> = U^/s, with 
s = r sin the cylindrical radius. The induction equation of the meridional magnetic potential a 
becomes: 

— ■ V(sArf)- Aifl = 0. (6) 
s 

In classical mean-field dynamo theory, one requires mechanisms to convert large-scale poloidal 
magnetic field into large-scale toroidal magnetic field, and vice versa. For rapidly rotating bod- 
ies, it is believed that azimuthal shear is an efficient way of achieving the former. It corresponds 
to the induction term of equation |5] and is therefore called the 'w-effect'. Our equation |6] for the 
meridional (poloidal) magnetic field only converts poloidal magnetic field (the imposed dipole) 
into poloidal magnetic field. One usually resorts to the contribution of small-scale fluctuations to 
produce an 'a-effect', which converts toroidal magnetic field into poloidal magnetic field. The 
objective of our study is to see how much can be explained by the large-scale flow alone, in order 
to measure the need for an a-effect in our experiment. 

4. Azimuthal fields: forward problem 

In this section, I want to express all the observables as functions of the azimuthal velocity 
field Uip, which I will invert for. t/^ projects onto the associated Legendre functions f' of order 
1 as: 

odd I 

where only odd degrees I are considered because t/^ is symmetric with respect to the equator. 

4.1. induction equation 

Similarly, b^ is written as: 



b^(r,0)^ ^ b„ir)Pl(cose), (8) 



where the sum is over even degrees n since we anticipate that is equatorially antisymmetric (I 
use n to denote even degrees and I for odd degrees). With this decomposition, the diffusion term 
of the induction equation (jSj becomes: 



1 dHrbn) n{n + 1) 



while the induction term is expressed as: 



b„ 



P,',(cos0), 



(9) 



r4 
+2 



odd I 

dUiir) 
dr 



de 

2U,(r) 



cos eP' (cos 61) 



(10) 



In order to keep all 6 dependence in the associated Legendre polynomials only, we make use of 
the following relationships: 



d0 



21 



cos P\ = 



21 



(11) 
(12) 



Injecting into equation (fTOt . the induction term becomes: 



{I - 3)U, + 2r 



dUi 



dUi 

-{I + 4)f/, + 2r— i 
dr 



dr 

{1+1)PU 



IP 



/+1 



(13) 



which confirms that only even degrees n contribute to the azimuthal magnetic field b^. Equating 
equation|9]and the opposite of equation[T3] we obtain for each even degree n a linear relationship 
between b„{r) and U„±\(r). Appendix B indicates the truncation degree to be considered for the 
various fields when l,„ax is the truncation degree of the expansion of the azimuthal velocity field 
in equation |7] 



4.2. electric field and potential 

In order to compare with the differences in electric potential measured at the surface of the 
sphere, we need to relate the electric potential V to the azimuthal velocity U^. The Eg component 
of the electric field E - - W can be decomposed as: 

Ee(r,e)^ ^ e„(r)pi(cos0), (14) 

even n 

which we integrate to get: 

y(r, 9) = V{r) - r e„(r)P°(cos 6). (15) 

even n 
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We compute the radial function V(r) by integrating Er along the z-axis. The Eg and E,- compo- 
nents of the electric field are derived from Ohm's law: 



E^j-UxB, (16) 

which yields: 

£«(r,0) = -J;K)-|X^'('-> 

odd I 

We get the e„{r) coefficients of equation [14] by identification, using the expansion of in b„, 
which were related to Ui in the previous paragraph. 



Z+ 1 1 I 



(17) 



4.3. boundary conditions 

We make use of the continuity of the magnetic field and of the tangential components of the 
electric field at the interface between the fluid and the solid shells. At these interfaces, thin veloc- 
ity boundary layers form in the fluid, which we will not be able to resolve from our Doppler pro- 



files. I therefore authorize azimuthal velocity discontinuities on the walls. Appendix C details 
how we deal with these discontinuities, as well as those in the electrical conductivity between the 
three shells. Note that electric currents are present in the solid shells and that the radial derivative 
of b^ is discontinuous. In the end, we obtain a relation between the magnetic field and its radial 
derivative (on the fluid side) at the interface r = r,- given by: 



db„ 
dr 



^-Pnbn(n)-Q„, (18) 

and similarly for r - rg. The expressions of !P„ and (3„ are given bv lC.lSl for r - r,. 



4.4. magnetic torque 

The magnetic torque is the torque due to the Lorentz force acting on the copper shell. We 
can compute it by integrating: 



M 



■ An 



Jf, Jo 



do r^ sin^ 6 TJr, 6), 



where is the azimuthal component of the Lorentz force per unit volume given by: 



— (sinO b^) + 2cos0 — (rba,) 
de ^ dr ^ 



(19) 



(20) 



Since we could project b^ over the P^{cos6), and because of the recursive relations given in 
equations [T2l the azimuthal component of the Lorentz force also projects onto the P^'(cos 6): 

r^(r,e)^Y,f'^''^Pi^'"''^^- (21) 

oddl 

Developing the expansion for b^p, we obtain: 



Mr) 



1 f 1 + 2 



\2l + 3 



{\-l)bM+2r 



dbi+i 

dr 
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/- 1 
'21- 1 



{I + 2)bi-i + 2r 



dbi-\ 
dr 



(22) 



Getting back to the expression of the magnetic torque (equation [T9]). we note that because 
of the orthogonaHty of the P'"(cos 9), the only degree of the expansion of the Lorentz force that 
contributes is / = 1, which only implies the degree 2 component of b^. It is then straightforward 
to see that the magnetic torque is directly proportional to the value of at the surface of the 
copper shell, and does not depend on its thickness. We simply get: 

Tm = -—biird- (23) 

Note that, more generally, the integral of the magnetic torque over the volume of a conductive 
shell reduces to a surface integral of the Maxwell tensor expressed as B,b^. 



5. Azimuthal inversion 



5.7. inversion set-up 

We now have all the expressions that relate our various measurements to the radial functions 
Ui{r) of the different degrees / of the expansion of the azimuthal velocity field. In order to carry 
out the inversion, I define a regular radial grid of Nr =137 points between r, = 74/210 and 
r,, - \, on which the Ui functions are discretized. The model vector M consists of Nr Umax + l)/2 
unknowns. 

The data vector D consists of a set of 5 azimuthal velocity profiles (with about 150 inde- 
pendent point measurements each) measured by ultrasound Doppler velocimetry, 16 azimuthal 
magnetic field measurements inside the fluid from 3 different latitudes, 4 electric potential dif- 
ferences at the surface, and the value of the torque on the inner sphere. The spatial data coverage 
is shown in figure IZXT] 



5.1.1. matrix formulation 

Radial derivatives are computed by finite difference. The diffusion operator in equation |9] 
is expressed in matrix form, taking into account the conditions at the boundaries expressed by 
equation [18] Under these conditions, the inverse problem is linear: all measured quantities are 
linearly related to the model vector, which consists of the f//(r^) values. Note that the contribution 
of the spinning inner sphere is treated separately, since there is no need to invert for this known 
(imposed) velocity. 

The best fitting model M is ob tained with the classical generalized least square inverse 



method dTarantola and Valettelll982l) : 



M = CppG^ [Cdd + GCppG'^y' D, (24) 

where G is the matrix from which the predicted data are computed by Dfn - GM. Cdd is the 
covariance matrix of the data, and Cpp the a priori covariance matrix of the model parameters, 
which I describe below. 



5.7.2. covariance matrices 

The covariance matrix of the data is taken diagonal. The diagonal terms are the square of 
the standard deviations presented in the data processing section |23] I use the a priori covariance 
matrix of the model to control the smoothness of the model in the radial direction, while the 
smoothness in the latitudinal direction is controlled by the truncation degree l,„ax of the Legendre 
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Figure 3: Azimuthal inversion: radial profiles of the modes of different harmonic degree. From left to right: radial 
functions of the azimuthal velocity f/;(r), the azimuthal magnetic field h„(r) and the orthoradial electric field e„(r). All 
fields are dimensionless as given in |Appendix A| The radius axis extends from fj to ?„. horizontal dotted lines indicate 
the fluid/solid interfaces at r, and r„. 



expansion. Within an /-block of the a priori covariance matrix of the parameters, the elements 
are taken as: 



Cnp(ra,rt) = 



To 


2 


I 


exp 



(25) 



where <5 is a smoothing parameter and tq/I is a damping parameter. 

In order to limit the velocity jump at the inner and outer spheres, where the Doppler profiles 
cannot resolve the boundary layers, we add fake profiles with zero velocity and a large error bar 

(T = 0.1). 



5.2. inversion results 

1 present the results of a simultaneous inversion of the velocity profiles, the induced magnetic 
field inside the fluid, the electric potential diff'erences at the surface, and the torque on the inner 
sphere. The data coverage is shown in figure 12.3.11 The Legendre expansion of the azimuthal 
velocity is carried up to Imaj, = 7 (/ odd), and the parameters of the a priori covariance matrix in 
expressionlZSl are chosen as 6 - 0.2 and tq - 0.01. 

5.2.7. radial profiles of the modes 

The radial profiles of the Legendre modes are shown in figure [3] for the coefficients of the 
azimuthal velocity f//, the azimuthal magnetic field b„ and the 0-component of the electric field 
e„. The I - 1 mode dominates the velocity field. The magnetic field modes show that the 
induction at the inner sphere boundary is strong, but that the shear inside the fluid produces large 
dipole in - 2) and octupole (« = 4) field as well. The electric field mode is almost entirely n = 2 
and does not drop to zero at the outer surface (r = rg). 



13 






Figure 4: Contour maps of the reconstructed angular velocity of the fluid (left), induced azimuthal magnetic field (center), 
and electric potential (light), from the inversion of DTS measurements. The amplitudes are given by the color bars. All 
fields are dimensionless as given in |Appendix A| 



5.2.2. maps of the fields 

Figure |4] shows contour maps (in the northern meridional quadrant) of the angular velocity 
of the fluid o), the azimuthal magnetic field b^, and the electric potential V. All quantities are 



dimensionless, as given in Appendix A 



The ma p of ang ular velocity illustrates some important properties of the flow already dis- 
cussed by iBrito et aP (1201 ih . Note the region near the equator of the inner sphere where the 
angular velocity of the fluid gets larger (1 < w < 1.3) than the a ngular velocity of the inner 
sphere. In that region, the fluid obeys Ferraro's law of isorotation (iFerraroi 119371) : the angular 
velocity is constant along magnetic field lines (white dashed lines). Further away from the inner 
sphere, the mean flow is geostrophic: the velocity does not vary along the rotation axis (vertical 
axis). Also note that the magnetic coupling between the inner copper shell and the fluid is very 
efficient and entrains the fluid at large angular velocities throughout the fluid volume. There is a 
region of marked weak velocity at the pole of the inner sphere, for which we have no explanation 
yet. 

The map of the induced azimuthal field illustrates that the field is mainly produced in two 
distinct regions: at the surface of the inner sphere, with alternating polarities, and in the large 
geostrophic shear zone. Note that the amplitudes range from -0.08 to 0.02 in the RmBo units we 
chose for and that Rm ^ 9.4 for / = 3 Hz. 

The electric potential map appears dominated by the geometry of the imposed dipole field. 



5.2.3. fits to the data 

Let's now see how well the reconstructed azimuthal velocity field fits the observations under 
the (jj-effect approximation. Figure |5] compares the predictions of our model with the velocity 
measurements along the ultrasound Doppler profiles, the azimuthal magnetic field inside the 
fluid, and the electric potential differences at the surface. In all cases, we compute the normalized 
misfit defined as: 



misfit — 




where dobs are the observed data with their error Tabs, and dpred the model prediction. 
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Figure 5: Azimuthal inversion: fits to tlie data. Ttiese plots compare tlie measured data, with their error bars, to the 
predictions of our inversion model for the azimuthal fields. We use the same colors and symbols as in figure imi deft), 
which gives the spatial distribution of the measurements. All data are dimensionless as given in [Appendix A| Left: 
angular velocity from Doppler velocimefi'y, as a function of the distance d along the ray of the ultrasound beam. Center: 
azimuthal magnetic field measured inside the fluid at different radii for 3 latitudes. Right: differences in electric potential 
between electrodes 10° degrees apart at the surface of the outer shell. 

The ultrasound Doppler profiles are well fit by our simple velocity model. The normalized 
misfit is 1.9. The fact that it is larger than one indicates that representation errors are present. 
Some of our assumptions (such as symmetry with respect to the equator) are partially violated. 
We observe for example that slightly different velocities are measured at points where Doppler 
profiles intersect in the {s,z) plane. 

The predicted induced azimuthal magnetic field under the w-effect approximation {i.e., taking 
into account the induction from the mean flow alone) falls in the right range, but deviates strongly 
from the observations at all latitudes for r < 0.6. The normalized misfit is 8. 

With a normalized misfit of 3.2, the differences in electric potential at the surface are not 
perfectly fit, but display the proper trend. 

The predicted magnetic torque is -0.27, slightly smaller than the observed value of -0.29 ± 
0.03. 

I find it quite remarkable that the simple velocity model I invert for provides a reasonable fit 
to observations as diverse as induced magnetic field, electric potential and torque, when plugged 
into the mean flow induction equation. By playing with the a priori covariance matrix, it is 
possible to improve the fit to the observations, at the expense of added complexity in the velocity 
field. However, the velocity field is very strongly constrained by our direct measurements of 
azimuthal velocity along criss-crossing profiles inside the fluid, and 1 have not been able to find 
models that could significantly improve the fit of the observed induced magnetic field inside 
the fluid. 
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6. Meridional fields: forward problem 



We now consider the poloidal (meridional) part of the velocity field. According to equation 
121 the poloidal velocity field is expressed as: 

Hp = V X ue^. (27) 

The direction of rotation of the meridional flow changes sign at the equator, which implies that 
the potential u projects on associated Legendre polynomials with even degree n: 

M(r,0)= 2 u„{r)Pl,(cos9). (28) 

even n 

We derive the expressions of the m, and uq components, which we use to relate to the velocity 
profiles measured by ultrasound Doppler velocimetry; 



Ur{r, 0) 



^ u„{r) 1 d I , X 

> sm0P'(cos0) 



u„{r) ^ du„(r) 
r dr 



(29) 
(30) 



6.1. induction equation 

Our next step is to derive a linear relation between the potential a of the magnetic field and 
the potential u of the flow from the induction equation |6] which I recall here; 



— ■ V(iArf)- Aja = 0. 
s 

The potential a projects onto the associated Legendre polynomials: 

a(r,6) = ai(r)Pj{cose). 

odd I 



(31) 



(32) 



where the sum is over odd degrees I since we anticipate that a is symmetric with respect to the 
equator With this decomposition, the diffusion term of the induction equationlSTIbecomes: 



Aia = ^ 



odd I 



1 d^irai) 
r dr^ 



-ai 



PjicosO), 



(33) 



while the induction term becomes: 



■ V(sAd) - —T (-Ur sin 6 + 2ug COS ff) . 

s 



(34) 



since the potential of the dipolar field is (sin 6')/r^. Injecting the expressions for Ur and ug 
(equations[30b. one gets the following expression: 

^■V(iA^) = --i Yj |««^(sin0P;,(cos0)) +2 L + r^J cos ^?Pi(cos 0)1 . (35) 
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Injecting the recursive relations [12] into the previous equation, we finally obtain: 



Up 1 1 

s n ^ 2n + 1 



(3 + «)m„ + 2r 

dr 



nPl 



(2 - n)Un + 2r— 
dr 



(36) 



which confirms that only odd degrees / contribute to the potential a of the meridional magnetic 
field bp. Equating equation[33]and equation|36] we obtain for each odd degree I a linear relation- 
ship between a;(r) and M/±i(r). [Appendix B recalls the truncation degree to be considered for 
the various fields when n„ax is the truncation degree of the expansion of the meridional velocity 
field potential u in equationl28] 



6.2. boundary conditions 

We impose non-penetration of the fluid at the interface with the solid shells. However, as for 
the azimuthal velocity, I allow for a discontinuity of the tangential velocity ue since we cannot 
resolve the very thin boundary layers. The continuity of the magnetic field implies that both a 
and its radial derivative are continuous. Only azimuthal currents contribute to the meridional 
magnetic field. Thus there is no electric current in the solid shells since both the azimuthal 
electric field and the U x B contribution to Ohm's law (equation \T6[ are null. It means that 
the meridional field matches a potential field on both sides of the fluid sheU, yielding a simple 
relation between a and its derivative: 

I a,(n) (37) 
-(/ + 1) a,(r,). (38) 



dai 
dr r, 
da I 
dr r„ 



7. Meridional inversion 



7.1. inversion set-up 

We now have all the expressions that relate our various measurements to the radial functions 
Un(r) of the different degrees n of the expansion of the meridional velocity potential. In order to 
carry out the inversion, I use the same regular radial grid of A^r = 137 points between r,- - 74/210 
and To = 1, on which the u„ functions are discretized. The model vector M consists of Nr n^^xl'^ 
unknowns. 

The data vector D consists of a set of 8 Doppler profiles (radial profiles of the radial velocity 
at 3 latitudes, one meridional profile, and meridional velocities retrieved from 4 off-meridional 
profiles by summing azimuthal profiles for opposite rotation direction), measurements of br and 
be inside the fluid at 3 latitudes, and measurements of br and bg at the surface from the equator 
up to 57° in latitude. The spatial data coverage is shown in figure IZXTI 

7.1.1. matrix formulation 

We follow the same approach as for the azimuthal inversion. Radial derivatives are computed 
by finite difference. The diffusion operator in equation|33]is expressed in matrix form, taking into 
account the conditions at the boundaries expressed by equation [38] Under these conditions, the 
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Figure 6: Meridional inversion: radial profiles of the modes of different harmonic degrees. From left to right: radial 
functions of the meridional velocity potential u„(r), and of the meridional magnetic potential a;(r). All fields are dimen- 
sionless as given in [Appendix A| The radius axis extends from r, to horizontal dotted lines indicate the fluid/solid 
interfaces at r, and . 



inverse problem is linear: all measured quantities are linearly related to the model vector, which 
consists of the u„{rii) values. The best fitting model M is obtained using the classical generalized 
least square inverse given by equationl24l 

Cdd is the covariance matrix of the data, and Cpp the a priori covariance matrix of the model 
parameters, which I describe below. 

7.1.2. covariance matrices 

The covariance matrix of the data is taken diagonal. The diagonal terms are the square of the 
standard deviations presented in the data processing section |273] The a priori covariance matrix 
of the model parameters has the same expression |25] as for the azimuthal inversion. In order to 
impose non-penetration of the fluid at the inner and outer shells, I force M„(r,) - Unir,,) - for 
all degrees n. 

7.2. inversion results 

I present the results of a simultaneous inversion of the velocity profiles and the induced 
magnetic field inside the fluid and at the surface. The data coverage is shown in figure IZTT] 
The Legendre expansion of the meridional velocity potential is carried up to n,„a,x = 8 (n even), 
and the parameters of the a priori covariance matrix in expre s sion |25] are chosen as (5 = 0.15 and 
TO =0.001. 

7.2.7. radial profiles of the modes 

The radial profiles of the Legendre modes are shown in figure |6] for the coefficients of the 
meridional potentials of velocity (m„) and magnetic field (a/) . All quantities are dimensionless. 



as given in Appendix A Note that the velocity modes show more structure in radius than the 



azimuthal velocity modes. 
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Figure 7: Stream lines of the meridional flow (left) and field lines of the induced meridional magnetic field (light) 
reconstructed from the inversion of DTS measurements. The flow is poleward beneath the outer shell. The values of the 
stream function range from -2.0 10"^ to 0.2 10"^ Those for the magnetic field from -1.0 10"^^ to 0.2 10"^ The blue 
sohd lines indicate anti-clock-wise circulation, and the pink dotted fines clok-wise circulation. The zero line is magenta. 



Note that the ind uced dipole field (/ - 1) dominates near the inner sphere, but vanishes at the 



surface as it should (ISpence et all 120061) . 



7.2.2. maps of the fields 

1 show the streamlines of the meridional flow, and the field lines of the induced poloidal 
magnetic field in figure Q 

The meridional circulation is poleward beneath the outer boundary, where th e lines are 



squeez ed, indicating large orthoradial velocities. Our solution is very similar to that of lBrito et al 



(1201 11) who inverted the velocity data alone, using a sine decomposition in radius. 

Not surprisingly, the induced magnetic field is largest near the inner sphere, where the im- 
posed dipole is strongest. The field lines are clockwise around the equator. Note that the induced 
field remains very small compared to the imposed dipole field, in agreement with the approxi- 
mation used in our forward problem (see section[3]). Indeed, I have plotted the field lines of the 
total field in figure |2T| (dashed lines), but they hardly depart from those of the dipole. 

7.2.3. fits to the data 

Let's now see how well the reconstructed meridional velocity field fits the observations. Fig- 
ure [8] compares the predictions of our model with the measured radial velocity profiles and the 
other meridional velocity profiles, all obtained by Doppler acoustics. The fit looks rather good, 
with an overall normalized misfit of 1.2. The more complex radial structure, as compared to the 
azimuthal case, seems required by the measured profiles. 

Figure |9] compares the observations to the predicted induced magnetic field, under the ap- 
proximation that it all results from induction by the mean meridional flow alone. The br and 
bg components inside the fluid are given in the first two panels. The model predicts reasonable 
amplitudes, which are in the range +0.01 in the RmBo units we chose for the induced field. 
Unfortunately, our data coverage is rather sparse, and does not extend deep into the fluid. The 
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Figure 8: Meridional inversion: fits to the velocity data. These plots compare the measured data, with their error bars, to 
the predictions of our inversion model for the meridional velocity field. We use the same colors as in figure [23T1 (right), 
which gives the spatial distribution of the measurements. All data are dimensionless as given in |Appendix A| Left: radial 
velocity from the radial Doppler velocimetry profiles at three latitudes, as a function of radius. Right: meridional velocity 
projected on the ultrasound beam as a function of the distance d along the beam for the non-radial profiles. 
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Figure 9: Meridional inversion: fits to tlie magnetic data. These plots compare the measured data, with their error bars, 
to the predictions of our inversion model for the induced meridional magnetic field. We use the same colors and symbols 
as in figure imi deft). which gives the spatial distribution of the measurements. All data are dimensionless as given in 
[Appendix A| Left: radial magnetic field measured inside the fluid at difl'erent radii for 3 latitudes. Center: orthoradial 
magnetic field measured inside the fluid at different radii for 3 latitudes. Right: radial and orthoradial components of the 
induced magnetic field measured at the surface of the outer shell. 



overall normalized misift is 2.4. The trends for br are compatible with the data, while there is a 
clear disagreement for bg. 

The right panel of figure|9]shows the predicted magnetic field at the surface of the outer shell. 
It fits rather nicely the data we measure, with a normalized misfit of 1 . 1 . 

As in the azimuthal case, I find it rather remarkable that the velocity model, which is strongly 
constrained by our Doppler profiles, predicts rather well our magnetic observations. It turns out 
that the magnetic field at the surface is quite sensitive to subtle changes in the velocity field. In 
the case I present, I over-emphasized the importance of the surface measurements to obtain a 
good fit. The slight modifications it brought to the velocity potential did not degrade the fit to the 
Doppler profiles. In contrast, I was not able to provide a better fit to the be field measured inside 
the fluid without introducing velocity features in contradiction with the velocity measurements. 
Remember however that all meridional observables are at least one order of magnitude less than 
their azimuthal counterpart, and that some biases cannot be totally excluded. 



8. Discussion 

In this article, I have exposed the formalism with which one can model the induced azimuthal 
magnetic field from the azimuthal velocity field (the w-effect), as well as other related measurable 
quantities (electric potentials and magnetic torque), taking into account realistic conductivity 
and velocity boundary conditions (section |4]i. I have written the inverse problem and inverted 
simultaneously ultrasound Doppler velocity profiles and magnetic observations from the DTS 
experiment, in order to recover the mean azimuthal flow (section|5]). I have carried out a similar 
procedure for the meridional flow and the magnetic field it induces (sections |6] and |7]l. 
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The results are rather remarkable, in that I obtain a full coherent solution for the mean ve- 
locity and magnetic fields inside the fluid (figures |4] and |7]i, which can be used to evaluate in 
a quantitative fashion the distribution of energies and dissipation, the force balance, and so on. 
In a sense, we have performed a numerical simulation of a magnetohydrodynamic flow at high 
Reynolds number (Re ^ 10*'), using the experiment as a Navier-Stokes solver, and working out 
the induction equation numerically. 

As an illustration, I report in table |2] the magnetic and kinetic energies of the azimuthal and 
meridional flows. They were obtained by integration of our solutions over the fluid shell (with 
the addition of the inner and outer solid shells for the magnetic field). All energies are scaled by 
pQ^{r*gf, yielding a value of order 1 for the kinetic energy of the azimuthal flow. Note that kinetic 
energies E"^ and e'^^'' are larger their magnetic counterparts e"^ and e"^'' . The magnetic energy of 
the imposed dipole field, whose energy, given by E^j''"''' = 51/Rm^, is in fact independent of the 
spin rate of the inner sphere. For the spin rate / = 3 Hz considered in this article (Rm ^ 9.4), the 
energy of the azimuthal flow E'l^ is comparable to that of the imposed magnetic field. 













sk 


51/Rm^ 


1.910-"* 


2.3 10-^ 


0.45 4 10-"^ 


~ 2 IQ-^ 


~ 4 10-^ 



Table 2: Energies of the imposed dipolar field, the azimuthal magnetic field and the meridional magnetic field, and of 
the azimuthal and meridional flows, all in units of pf2^(r*)' = prj ^r'Rm^. The table als o lists estimates of the energies 
of the magnetic and velocity fluctuations Sm and ex deduced from lFigueroa et al] )2012h in the same units. 



The mean flow and magnetic field in the DTS experiment have been examined bv lBrito et al 



( 1201 Ih . who show that at first order both the azimuthal and meridional velocities and the induced 
azimuthal magnetic field increase linearly with the rotation frequency / of the inner sphere (or 
equivalently with the magnetic Reynolds number Rm) for / up to 30 Hz. The induced meridional 
magnetic field evolves in a more complex fashion, but 1 will assume here that it remains linear 
in Rm as well. As a consequence, the fields we have inverted for at / = 3 Hz give a rough first 
order image of the fields at higher rotation frequency. 

Mean-field theory indicates that correlated fluctuations of the velocity and magnetic fields 
can contribute to the axisymm etric mean fields we are d ealing with. Fluctuati o ns in the DTS- 



experi ment are described in ,Schmitt et al.l (120081 120131) : iFigueroa et alJ (l2012h . iFigueroa et al 



(120 12h show that velocity and magnetic fluctuations also increase linearly with Rm at first order, 
and give an estimate for the energy of the fluctuations, which is listed in table|2]for comparison. 

1 will use my inversion results and these various estimates to evaluate the contribution of the 
terms that have been left out of the induction equation under my simplifying assumptions. The 
complete induction equation for the axisymmetric steady state, which replaces equations |5] and 
|6l can be written as: 



A,/7, 



(39) 



-sBd ■ -I- Rm ^-sbp ■ Vcj + sUp-V (^) - [v x <m x ft)]^ 

All symbols have the same meaning and adimensionalization as in the rest of this article, yielding 
the Rm-dependence of the additional terms. Let us now check that these additional terms are 
negligible indeed at the rotation rate we considered (/ = 3 Hz, which yields Rm = 9.4), and 
evaluate their influence at the highest rotation rate we achieve (/ = 30 Hz, yielding Rm - 94). 
Before doing so, we examine the various additional terms inside the braces of equation[39l 
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Figure 10: (left and center) maps of the contributions of the additional induction terms from equations l39l to the azimuthal 
magnetic field and to the meridional field bp. The values of the magnetic stream function for bp range from -1.2 10"* 
to 0.5 10"'. (right) tentative map of the product of average fluctuations in velocity and magnetic field (u){b). 



Compared with equation|5]used for the inversion, the equation for contains three additional 
terms. The first one corresponds to the w-efFect applied to the induced meridional magnetic 
field. The second one is the advection of the azimuthal field by the meridional flow. The third 
one describes the mean-field axisymmetric induction caused by the correlated non-axisymmetric 
fluctuations of velocity u and magnetic field b. We can compute the first two coupling terms by 
combining the results of our azimuthal and meridional inversions. The products are computed 
in the physical space, and projected onto the P,' Legendre polynomials (see Appendix B for 
the parity and bounds of the projection). Solving the induction equation for these two coupling 
terms, I get the map of the azimuthal magnetic field they produce (figure [TOt. Its amplitude 
amounts to about 10""* Rm, yielding a contribution of less than l%tob^ for Rm = 9.4 (compare 
with figure |4]i. 

Only one coupling induction term appears in the equation for the meridional magnetic poten- 
tial a. It corresponds to the advection of the induced meridional field by the meridional flow. I 
compute it and solve the induction equation to get the field lines of its contribution to the merid- 
ional magnetic field, shown in figure [10] Its amplitude is about 10"'' Rm, which modifies the 
inverted meridional magnetic field of figure Qby about 1% for Rm = 9.4. 

I now turn to the last term of both equations, which is the most interesting one because 
it describes the contribution of non-axisymmetric fluctuations to the mean axisymmetric mag- 
netic field. This effect is supposed to play a major role in the generation of large-scale mag- 
netic fields by dynamo action. There exists a vast literature on the axisymmetric mean-field 
dynamo equations, which have proven ver y useful in de ciphering the magnetic behaviour of 
the Sun and stars. The excellent review bv ICharbonneaul (12005) recalls that the electromotive 
force & - {u xb) term is often expressed as a truncated series expansion of the large-scale 
magnetic field: & - a : {B) + ^ . Y x {B}, where a and )8 are two tensors, which depend 
upon the turbulent characteristics of the flow. The c-effect is crucial for transferring energy 
from the azimuthal magnetic field to the meridional field, while the y6-effect can be seen as en- 
hancing the magnetic diflFusivity. Although some theories can provide e xpressions for these two 
effects in idealized situations such as homogeneous isotropic turbulence (ISteenbeck et al. ■ Il966l) 



or magnetohydrodynamic turbulence ( Pouquet et al.[[l9 76). or phase fluctuations of simple flows 
dPetrelis and Fauve , l2006h . it remains a challenge to estimate their contribution in actual experi- 
ments and in celestial objects. Recent experiments have shed some light on this issue: 
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The torus experiment in Perm consists in spinning at rotation rates up to 45 Hz a torus filled 
with liquid sodium. The torus is then stop ped abruptly, generating a highly turbulent decaying 
flow. Using this device, iFrick et al. I (l2010h have measured the electromagnetic response of the 
torus channel to an alternating magnetic field, thereby measuring the eff'ective global magnetic 
diff'usivity. They have observed an enhanced diffusivity during the turbulent decay, which they at- 
tribute to a yS-effect. The magnetic diff'usivity was found to be up to 30% larger than the molecular 
diffusivity, for the highest magnet ic Reynolds number Rm ^ 30. Local velocities were measured 
in the same experiments, enabling lNoskov et al.l(l2012b to compare turbulent magnetic diffusivity 
with turbulent viscosity, and to propose a scaling for turbulent diffusivity as: p ~ Rm^,„5, where 
RiHrmj = Urmsro/t] is a magnetic Reynolds number based on the amplitude of the velocity fluctu- 
ations Urms, and To is the channel radius. The exponent k is about 2 for Rm„„j < 1, decreasing to 
about 1 .3 above that. A maximum value of Rm,„,5 ^ 1 .5 was achieved in this experiment. 

The Madison liquid sodium experiment has two counter-rotating impellers facing each other 
in a sphere filled with 520 liters of liquid sodium . With a maximum r otation rate of 30 Hz, mag- 
netic Reynolds numbers up to 160 are achieved. ISpence et al. I (l2006l) have reported a turbulent- 
induced large-scale magnetic field in this set-up. Fluctuations induce some a-effect, which de- 
creases the in duction produced by the mean large-scale flow, thereby preventing dynamo onset. 
More recentlv. lRahbarnia et al.l (120 12h have installed a probe that measures the velocity and mag- 
netic field vectors simultaneously at a single position. This gives access to components aligned 
with the mean magnetic fi eld (akin to an g-effec t) and to those aligned with the mean electric 
current (akin to a yS-effect). iRahbarnia et al. ( 2012 ) find that the latter dominates in the current set- 
up that includes an equatorial baffle, which reduces the amplitude of the largest-scale turbulent 
eddies. For their highest magnetic Reyn olds number Rm"^^ - 160, they measure an excess mag- 
netic diffusivity of about yS - 30%. As iNoskov et al. I (l2012h they introduce a fluctuation-based 
magnetic Reynolds number: Rm',,,, - Urmsth, which differs because the typical dimension ( is 
the correlation length of the fluctuations, which they find to be at least one order of magnitude 
smaller than the integral scale. The maximum value they reach is Rm*^,, ^ 1 .2, corresponding to 
Rm„„j - 12. The relation [5 ^ Rm*„,j looks compatible with their data. 

Getting back to DTS and to equations|39] I note that, unlike the coupling terms, the induction 
due to non-axisymmetric fluctuations cannot be computed from my inversion results. In fact, we 
have no access to the {u x h) term. Howeve r, we can estima t e a rea sonable upper bound from 
the product {u)(h). In the example shown bv IRahbarnia et al. ( 2012 ). the average electromotive 
force amounts to about 60% of this product. Figure [lO] displays a tentative map of {u){b) in 
the DTS experiment. Unlike the other maps in this article, it does not result from a rigorous 
inversion. I have assumed a decomposition on (with even n) of the average velocity and 
magnetic fluctuations. I thus constructed a map of the velocity fluctuations from the inversion 
of the root-mean square fluctuations measured along the same Doppler profiles as for the mean- 
flow inversions. Similarly, a map of the magnetic fluctuations was built by inversion of the rms 
fluctuations of the magnetic field in the sleeve at the three different latitudes. The map shown in 
figure[TO]is the product of these two maps. 

The lack of directional and time-correlation information precludes using this map to solve 
for the induced magnetic field directly, but we can compare the amplitude of the estimate of this 
induction term with that of the other additional terms discussed above. The amplitude of {u){b) 
reaches about 10""* Rm, while the amplitude of the additional coupling terms is about 2 10"^ Rm 
in the equation and 5 10 Rm in the equation for a. Therefore, non-axisymmetric fluctuations 
can contribute at best to a few percents of the total induction at our highest magnetic Reynolds 
number (Rm ^ 100), and should always remain smaller than the terms that couple the meridional 
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and azimuthal flows and magnetic fields. 

It may look surprising that fluc tuations have s o little effect in our case, while we reach mag- 
netic Reynolds number larger than lNoskov et alj (120121) . We know that velocity fluctuati ons are 



hampered by the presence of the strong imposed magnetic field (Figuer oa et al.L 120121) . Nev- 
ertheless, using our estimate of ck in table |2] we get Rm,,,,, ^ 4.5 10"^ Rm, yielding a value 
R mrmi - 4 at our hig hest rotation rate. This is at least twice larger th an the maximum value 
of iNoskov et al.l (120121) . but comparable to the lowest value reported bv lRahbarnia et al.l (120 12h 
(Rm™, ^ 3.4 for Rm"'' = 60), for which the yS-effect is only 2%. 

However, one should keep in mind that I have focused on the inversion of data obtained at 
low magnetic Reynolds nu mber (Rm = 9.4), a nd assumed that the results can be extrapolated to 
our maximum Rm because iBrito et al.l (1201 lb observe that velocity and induced magnetic field 
are proportional to Rm at first order. But they also observe some important diff'erences when Rm 
increases and causes the relative strength of the imposed dipolar field to weaken: the zone of 
super-rotation shrinks and the intensity of the induced field increases more rapidly with Rm. We 
will have to investigate the role of these modifications. It remains that I was not able to perfectly 
fit all data at Rm = 9.4. In particular, my model does not predict high enough values for b^p, while 
it overestimates the electric potential differences measured at the surface. It could due to the fact 
that the experiment deviates somewhat from the assumptions used: the imposed magnetic field is 
not perfectly dipolar, the inner sphere is held by shafts, the equatorial symmetry could be slightly 
violated, etc. 



Nataf and Gagnierd (120081) note that fluctuations get even smaller when both the Coriolis 



and the Lorentz forces are dominant, as evidenced from measurements in the DTS experiment 
when the outer sphere spins. It would be interesting to investigate the role of fluctuations when 
these constraints are still present but at a much larger magnetic Reynolds number, such as can be 
achieved in the 3m-diameter BigSister experiment of Dan Lathrop and his group at the University 
of Maryland. 
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Appendix A. Adimensionalization 



quantity scale 

time n ' 

length r* 

velocity r^Q. - ^Rm 

imposed magnetic field Bg 

induced magnetic field Rm Bg 

electric field Rm rjBoIrl 

electric potential Rm tjBq 

electric current density Rm rjcrBo/rl 

magnetic torque Rm (rD^Bl/fio 
energy pQ^(r*)^ - prir*Rm^ 

Table A. 3: Adimensionalization used in this article, n is the imposed angular velocity of the inner sphere, and r* is 
the dimensional inner radius of the outer shell. is the intensity of the imposed dipolar magnetic field at the equator 
for r = Tq. Rm is the magnetic Reynolds number defined as Rm = (r*)^ i^/l- The magnetic diffusivity, the electric 
conductivity, and the density of the liquid are noted cr and p, respectively. 



Appendix B. Parity and truncation 



(r,^) field degree parity minimum degree maximum degree 







odd 


1 






by, 


even 


2 


^max "I" 1 




Eg 


even 


2 


^max "I" 1 






odd 


1 


^max "1" 2 




u 


even 


2 


^max 




a 


odd 


1 


^max 1 


coupling 


terms for 


even 


2 


^max ~^ ^max 1 


coupling 


I terms for a 


odd 


1 


2 * f^max ~ 1 



Table B.4: Parity and truncation degree of the various (r,0) fields analyzed in this article. Azimuthal inversion: U^, 
Eg and T^; meridional inversion: u and a; coupling terms in equation |39]f or and a. 



Appendix C. Dealing with discontinuities in velocity and conductivity 

In this appendix, I give the analytic expressions of the magnetic and electric fields in the 
solid spherical shells that contain the liquid sodium. I take into account the different electric 
conductivities and the contribution of a discontinuity of the tangential velocity at the interface 
between the solid and the fluid. I take advantage of the fact that we only consider the steady-state 
solution. 
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Appendix C.l. magnetic field 

In the solid inner and outer shells, equation ^ reduces to Ai/?^ - 0, which implies that the 
bn{r) functions of the decomposition ([8]) satisfy: 



1 d^(rb„) n(n + 1) 
r dr^ 



b„ = 0, 



(C.l) 



whose solutions are written as: 



'b„{r) = c^r"' +c„r 



(n+l) 



(C.2) 



where the * left superscript refers to the solid shell. Since electric currents cannot circulate 
outside the shells, b^p - Q ai r - f-j and r - r,,. Adding the continuity of b^ across the boundaries, 
we obtain the following relations between the ^ and " coefficients of equation (IC.2l l: 



c r: + c r- 



= 

= b„(n). 



Solving this linear system, we obtain: 



-(r7(«+'>/D„) b„(ri) 



c- = (ri/D„)b„(n), 



where the determinant D„ is: 



-(«+i) 



(C.3) 

(C.4) 
(C.5) 

(C.6) 



The continuity of the 9 component of the electric field yields the final relation. On the two 
sides of the solid/fluid interface, we have: 



(C.7) 



where the ' left superscript refers to the liquid. The continuity of Eg thus yields the following 
relation between the jumps in velocity and conductivity, and that in the derivative of the magnetic 
field at r, (with our adimensionalization given in Appendix A i: 

d(rb^) ^»cr d(r^"b^) 




rdr 



CUf 



rdr 



-[U^{ri,e)- C"f/^(r,,0)]B,(r,-,0), 



(C.8) 



where ^"cr and ''"cr are the electrical conductivities of liquid sodium and of the copper inner shell, 
respectively. The azimuthal velocity of the inner sphere is "^"[/^(r,, 6*) = r, sinfl = -r, Pj(cos ^?). 



We now evaluate the derivative term of equations IC.Sl in the solid part, using expression lC.2 



rdr 



= Yj Cnb„(n)Pl{cOS0), 



with: 



-(n+2) 

■nr. r" 



(C.9) 



(C.IO) 
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Injecting equations IC. 9 1 in |C8] and expressing the radial derivative of rb^ in the fluid, we get: 



C„\bn{ri) + 



db„ 
dr 



P,',(cos 6)^- [U^iri, 6) - n sin o] BAn, 6). (C.l 1) 



Making use of the Legendre projection ( fTTl i for the ux B term, we finaUy obtain: 



1 



Na 



(T 



i odd I 



cr 
/+ 1 



db„ 



21 +\ 



-P , + 



I 



21+1 



(C.12) 



where dn is the Kronecker symbol. By identification of the of identical degree, we deduce a 
relationship between the magnetic field and its derivative at the interfaces for each even degree 

2<n< l,„ax + 1: 



db„ 
~d7 



-P„b„(n)-Qn, 



with: 



Qn 



--—Cn 

r; ''"cr 

2 / n + 2 n-\ 



(C.13) 

(C.14) 
(C.15) 



This condition is similar (albeit more complex) to the one encountered for the magnetic 
potential induced by the meridional circulation. The boundary condition at the interface between 
the fluid and the outer shell is exactly the same, replacing r, by and r, by fg, '""cr by cr (the 
electric conductivity of the stainless steel outer shell), and omitting the 6„2 term in Q„ since the 
outer sphere is at rest. 



Appendix C.2. electric field 

Using the expression IC.7I of the 9 component of the electric field in the solid shell, and 
injecting the decomposition lC.2l of the bn functions, we obtain an analytical expression for the e„ 
functions in the inner copper shell: 

c«e„{r) = ^ [(« + l)r:("^^V«-' + n f] r-^-^)] M^ll + ^ ^^^^^ 

A similar relation is easily derived for the outer stainless steel shell. 
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